ProNet DB: a proteome-wise database for protein surface property representations and RNA-binding profiles.

The rapid growth in the number of experimental and predicted protein structures and more complicated protein structures poses a significant challenge for computational biology in leveraging structural information and accurate representation of protein surface properties. Recently, AlphaFold2 released the comprehensive proteomes of various species, and protein surface property representation plays a crucial role in protein-molecule interaction predictions, including those involving proteins, nucleic acids and compounds. Here, we proposed the first extensive database, namely ProNet DB, that integrates multiple protein surface representations and RNA-binding landscape for 326 175 protein structures. This collection encompasses the 16 model organism proteomes from the AlphaFold Protein Structure Database and experimentally validated structures from the Protein Data Bank. For each protein, ProNet DB provides access to the original protein structures along with the detailed surface property representations encompassing hydrophobicity, charge distribution and hydrogen bonding potential as well as interactive features such as the interacting face and RNA-binding sites and preferences. To facilitate an intuitive interpretation of these properties and the RNA-binding landscape, ProNet DB incorporates visualization tools like Mol* and an Online 3D Viewer, allowing for the direct observation and analysis of these representations on protein surfaces. The availability of pre-computed features enables instantaneous access for users, significantly advancing computational biology research in areas such as molecular mechanism elucidation, geometry-based drug discovery and the development of novel therapeutic approaches. Database URL:  https://proj.cse.cuhk.edu.hk/aihlab/pronet/.

Monitoring Real-Time Temperature Dynamics of a Short RNA Hairpin Using Förster Resonance Energy Transfer and Circular Dichroism.

RNA molecules play crucial roles in gene expression regulation and cellular signaling, and these functions are governed by the formation of RNA secondary and tertiary structures. These structures are highly dynamic and subject to rapid changes in response to environmental cues, temperature in particular. Thermosensitive RNA secondary structures have been harnessed by multiple organisms to survey their temperature environment and to adjust gene expression accordingly. It is thus highly desirable to observe RNA structural changes in real time over a range of temperatures. Multiple approaches have been developed to study structural dynamics, but many of these require extensive processing of the RNA, large amounts of RNA input, and/or cannot be applied under physiological conditions. Here, we describe the use of a dually fluorescently labeled RNA oligonucleotide (containing a predicted hairpin structure) to monitor subtle RNA structural dynamics in vitro by Förster resonance energy transfer (FRET) and circular dichroism (CD) spectroscopy. These approaches can be employed under physiologically relevant conditions over a range of temperatures and with RNA concentrations as low as 200 nM; they enable us to observe RNA structural dynamics in real time and to correlate these dynamics with changes in biological processes such as translation.

Global and single-nucleotide resolution detection of 7-methylguanosine in RNA.

RNA modifications, including N-7-methylguanosine (m7G), are pivotal in governing RNA stability and gene expression regulation. The accurate detection of internal m7G modifications is of paramount significance, given recent associations between altered m7G deposition and elevated expression of the methyltransferase METTL1 in various human cancers. The development of robust m7G detection techniques has posed a significant challenge in the field of epitranscriptomics. In this study, we introduce two methodologies for the global and accurate identification of m7G modifications in human RNA. We introduce borohydride reduction sequencing (Bo-Seq), which provides base resolution mapping of m7G modifications. Bo-Seq achieves exceptional performance through the optimization of RNA depurination and scission, involving the strategic use of high concentrations of NaBH4, neutral pH and the addition of 7-methylguanosine monophosphate (m7GMP) during the reducing reaction. Notably, compared to NaBH4-based methods, Bo-Seq enhances the m7G detection performance, and simplifies the detection process, eliminating the necessity for intricate chemical steps and reducing the protocol duration. In addition, we present an antibody-based approach, which enables the assessment of m7G relative levels across RNA molecules and biological samples, however it should be used with caution due to limitations associated with variations in antibody quality between batches. In summary, our novel approaches address the pressing need for reliable and accessible methods to detect RNA m7G methylation in human cells. These advancements hold the potential to catalyse future investigations in the critical field of epitranscriptomics, shedding light on the complex regulatory roles of m7G in gene expression and its implications in cancer biology.

Expedient production of site specifically nucleobase-labelled or hypermodified RNA with engineered thermophilic DNA polymerases.

Innovative approaches to controlled nucleobase-modified RNA synthesis are urgently needed to support RNA biology exploration and to synthesize potential RNA therapeutics. Here we present a strategy for enzymatic construction of nucleobase-modified RNA based on primer-dependent engineered thermophilic DNA polymerases - SFM4-3 and TGK. We demonstrate introduction of one or several different base-modified nucleotides in one strand including hypermodified RNA containing all four modified nucleotides bearing four different substituents, as well as strategy for primer segment removal. We also show facile site-specific or segmented introduction of fluorophores or other functional groups at defined positions in variety of RNA molecules, including structured or long mRNA. Intriguing translation efficacy of single-site modified mRNAs underscores the necessity to study isolated modifications placed at designer positions to disentangle their biological effects and enable development of improved mRNA therapeutics. Our toolbox paves the way for more precise dissecting RNA structures and functions, as well as for construction of diverse types of base-functionalized RNA for therapeutic applications and diagnostics.

Catalytic Metal Ion-Substrate Coordination during Nonenzymatic RNA Primer Extension.

Nonenzymatic template-directed RNA copying requires catalysis by divalent metal ions. The primer extension reaction involves the attack of the primer 3'-hydroxyl on the adjacent phosphate of a 5'-5'-imidazolium-bridged dinucleotide substrate. However, the nature of the interaction of the catalytic metal ion with the reaction center remains unclear. To explore the coordination of the catalytic metal ion with the imidazolium-bridged dinucleotide substrate, we examined catalysis by oxophilic and thiophilic metal ions with both diastereomers of phosphorothioate-modified substrates. We show that Mg2+ and Cd2+ exhibit opposite preferences for the two phosphorothioate substrate diastereomers, indicating a stereospecific interaction of the divalent cation with one of the nonbridging phosphorus substituents. High-resolution X-ray crystal structures of the products of primer extension with phosphorothioate substrates reveal the absolute stereochemistry of this interaction and indicate that catalysis by Mg2+ involves inner-sphere coordination with the nonbridging phosphate oxygen in the pro-SP position, while thiophilic cadmium ions interact with sulfur in the same position, as in one of the two phosphorothioate substrates. These results collectively suggest that during nonenzymatic RNA primer extension with a 5'-5'-imidazolium-bridged dinucleotide substrate the interaction of the catalytic Mg2+ ion with the pro-SP oxygen of the reactive phosphate plays a crucial role in the metal-catalyzed SN2(P) reaction.

ABC2A: A Straightforward and Fast Method for the Accurate Backmapping of RNA Coarse-Grained Models to All-Atom Structures.

RNAs play crucial roles in various essential biological functions, including catalysis and gene regulation. Despite the widespread use of coarse-grained (CG) models/simulations to study RNA 3D structures and dynamics, their direct application is challenging due to the lack of atomic detail. Therefore, the reconstruction of full atomic structures is desirable. In this study, we introduced a straightforward method called ABC2A for reconstructing all-atom structures from RNA CG models. ABC2A utilizes diverse nucleotide fragments from known structures to assemble full atomic structures based on the CG atoms. The diversification of assembly fragments beyond standard A-form ones, commonly used in other programs, combined with a highly simplified structure refinement process, ensures that ABC2A achieves both high accuracy and rapid speed. Tests on a recent large dataset of 361 RNA experimental structures (30-692 nt) indicate that ABC2A can reconstruct full atomic structures from three-bead CG models with a mean RMSD of ~0.34 Å from experimental structures and an average runtime of ~0.5 s (maximum runtime < 2.5 s). Compared to the state-of-the-art Arena, ABC2A achieves a ~25% improvement in accuracy and is five times faster in speed.

Architecture of an RNA polymerase ribozyme illuminates the RNA World.

'Ribo-organisms' of the primordial RNA World would have needed ribozymes that catalyze RNA replication. McRae, Wan, Kristoffersen et al. recently revealed how these RNA replicases might have functioned by solving the structure of an artificial polymerase ribozyme. This work illustrates how complex RNA structures evolve, with implications for the origins of life.

Spotlight on G-Quadruplexes: From Structure and Modulation to Physiological and Pathological Roles.

G-quadruplexes or G4s are non-canonical secondary structures of nucleic acids characterized by guanines arranged in stacked tetraplex arrays. Decades of research into these peculiar assemblies of DNA and RNA, fueled by the development and optimization of a vast array of techniques and assays, has resulted in a large amount of information regarding their structure, stability, localization, and biological significance in native systems. A plethora of articles have reported the roles of G-quadruplexes in multiple pathways across several species, ranging from gene expression regulation to RNA biogenesis and trafficking, DNA replication, and genome maintenance. Crucially, a large amount of experimental evidence has highlighted the roles of G-quadruplexes in cancer biology and other pathologies, pointing at these structurally unique guanine assemblies as amenable drug targets. Given the rapid expansion of this field of research, this review aims at summarizing all the relevant aspects of G-quadruplex biology by combining and discussing results from seminal works as well as more recent and cutting-edge experimental evidence. Additionally, the most common methodologies used to study G4s are presented to aid the reader in critically interpreting and integrating experimental data.

Coarse-grained modeling of DNA-RNA hybrids.

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

Chemoenzymatic Installation of Site-Specific Chemical Groups on DNA Enhances the Catalytic Activity.

Functional DNAs are valuable molecular tools in chemical biology and analytical chemistry but suffer from low activities due to their limited chemical functionalities. Here, we present a chemoenzymatic method for site-specific installation of diverse functional groups on DNA, and showcase the application of this method to enhance the catalytic activity of a DNA catalyst. Through chemoenzymatic introduction of distinct chemical groups, such as hydroxyl, carboxyl, and benzyl, at specific positions, we achieve significant enhancements in the catalytic activity of the RNA-cleaving deoxyribozyme 10-23. A single carboxyl modification results in a 100-fold increase, while dual modifications (carboxyl and benzyl) yield an approximately 700-fold increase in activity when an RNA cleavage reaction is catalyzed on a DNA-RNA chimeric substrate. The resulting dually modified DNA catalyst, CaBn, exhibits a kobs of 3.76 min-1 in the presence of 1 mM Mg2+ and can be employed for fluorescent imaging of intracellular magnesium ions. Molecular dynamics simulations reveal the superior capability of CaBn to recruit magnesium ions to metal-ion-binding site 2 and adopt a catalytically competent conformation. Our work provides a broadly accessible strategy for DNA functionalization with diverse chemical modifications, and CaBn offers a highly active DNA catalyst with immense potential in chemistry and biotechnology.

Enzymatic Assembly of Small Synthetic Genes with Repetitive Elements.

Gene synthesis efficiency has greatly improved in recent years but is limited when it comes to repetitive sequences, which results in synthesis failure or delays by DNA synthesis vendors. This represents a major obstacle for the development of synthetic biology since repetitive elements are increasingly being used in the design of genetic circuits and design of biomolecular nanostructures. Here, we describe a method for the assembly of small synthetic genes with repetitive elements: First, a gene of interest is split in silico into small synthons of up to 80 base pairs flanked by Golden-Gate-compatible overhangs. Then, synthons are made by oligo extension and finally assembled into a synthetic gene by Golden Gate Assembly. We demonstrate the method by constructing eight challenging genes with repetitive elements, e.g., multiple repeats of RNA aptamers and RNA origami scaffolds with multiple identical aptamers. The genes range in size from 133 to 456 base pairs and are assembled with fidelities of up to 87.5%. The method was developed to facilitate our own specific research but may be of general use for constructing challenging and repetitive genes and, thus, a valuable addition to the molecular cloning toolbox.

Atomic Mutagenesis of N6-Methyladenosine Reveals Distinct Recognition Modes by Human m6A Reader and Eraser Proteins.

N6-methyladenosine (m6A) is an important modified nucleoside in cellular RNA associated with multiple cellular processes and is implicated in diseases. The enzymes associated with the dynamic installation and removal of m6A are heavily investigated targets for drug research, which requires detailed knowledge of the recognition modes of m6A by proteins. Here, we use atomic mutagenesis of m6A to systematically investigate the mechanisms of the two human m6A demethylase enzymes FTO and ALKBH5 and the binding modes of YTH reader proteins YTHDF2/DC1/DC2. Atomic mutagenesis refers to atom-specific changes that are introduced by chemical synthesis, such as the replacement of nitrogen by carbon atoms. Synthetic RNA oligonucleotides containing site-specifically incorporated 1-deaza-, 3-deaza-, and 7-deaza-m6A nucleosides were prepared by solid-phase synthesis and their RNA binding and demethylation by recombinant proteins were evaluated. We found distinct differences in substrate recognition and transformation and revealed structural preferences for the enzymatic activity. The deaza m6A analogues introduced in this work will be useful probes for other proteins in m6A research.

Efficient Prediction Model of mRNA End-to-End Distance and Conformation: Three-Dimensional RNA Illustration Program (TRIP).

The secondary and tertiary structures of RNA play a vital role in the regulation of biological reactions. These structures have been experimentally studied through in vivo and in vitro analyses, and in silico models have become increasingly accurate in predicting them. Recent technologies have diversified RNA structure predictions, from the earliest thermodynamic and molecular dynamic-based RNA structure predictions to deep learning-based conformation predictions in the past decade. While most research on RNA structure prediction has focused on short non-coding RNAs, there has been limited research on predicting the conformation of longer mRNAs. Our study introduces a computer simulation model called the Three-dimensional RNA Illustration Program (TRIP). TRIP is based on single-chain models and angle restriction of each bead component from previously reported single-molecule fluorescence in situ hybridization (smFISH) experiments. TRIP is a fast and efficient application that only requires up to three inputs to acquire outputs. It can also provide a rough visualization of the 3D conformation of RNA, making it a valuable tool for predicting RNA end-to-end distance.

Influence of ion and hydration atmospheres on RNA structure and dynamics: insights from advanced theoretical and computational methods.

RNA, a highly charged biopolymer composed of negatively charged phosphate groups, defies electrostatic repulsion to adopt well-defined, compact structures. Hence, the presence of positively charged metal ions is crucial not only for RNA's charge neutralization, but they also coherently decorate the ion atmosphere of RNA to stabilize its compact fold. This feature article elucidates various modes of close RNA-ion interactions, with a special emphasis on Mg2+ as an outer-sphere and inner-sphere ion. Through examples, we highlight how inner-sphere chelated Mg2+ stabilizes RNA pseudoknots, while outer-sphere ions can also exert long-range electrostatic interactions, inducing groove narrowing, coaxial helical stacking, and RNA ring formation. In addition to investigating the RNA's ion environment, we note that the RNA's hydration environment is relatively underexplored. Our study delves into its profound interplay with the structural dynamics of RNA, employing state-of-the-art atomistic simulation techniques. Through examples, we illustrate how specific ions and water molecules are associated with RNA functions, leveraging atomistic simulations to identify preferential ion binding and hydration sites. However, understanding their impact(s) on the RNA structure remains challenging due to the involvement of large length and long time scales associated with RNA's dynamic nature. Nevertheless, our contributions and recent advances in coarse-grained simulation techniques offer insights into large-scale structural changes dynamically linked to the RNA ion atmosphere. In this connection, we also review how different cutting-edge computational simulation methods provide a microscopic lens into the influence of ions and hydration on RNA structure and dynamics, elucidating distinct ion atmospheric components and specific hydration layers and their individual and collective impacts.

HIPPO: HIstogram-based Pseudo-POtential for scoring protein-ssRNA fragment-based docking poses.

BACKGROUND: The RNA-Recognition motif (RRM) is a protein domain that binds single-stranded RNA (ssRNA) and is present in as much as 2% of the human genome. Despite this important role in biology, RRM-ssRNA interactions are very challenging to study on the structural level because of the remarkable flexibility of ssRNA. In the absence of atomic-level experimental data, the only method able to predict the 3D structure of protein-ssRNA complexes with any degree of accuracy is ssRNA'TTRACT, an ssRNA fragment-based docking approach using ATTRACT. However, since ATTRACT parameters are not ssRNA-specific and were determined in 2010, there is substantial opportunity for enhancement. RESULTS: Here we present HIPPO, a composite RRM-ssRNA scoring potential derived analytically from contact frequencies in near-native versus non-native docking models. HIPPO consists of a consensus of four distinct potentials, each extracted from a distinct reference pool of protein-trinucleotide docking decoys. To score a docking pose with one potential, for each pair of RNA-protein coarse-grained bead types, each contact is awarded or penalised according to the relative frequencies of this contact distance range among the correct and incorrect poses of the reference pool. Validated on a fragment-based docking benchmark of 57 experimentally solved RRM-ssRNA complexes, HIPPO achieved a threefold or higher enrichment for half of the fragments, versus only a quarter with the ATTRACT scoring function. In particular, HIPPO drastically improved the chance of very high enrichment (12-fold or higher), a scenario where the incremental modelling of entire ssRNA chains from fragments becomes viable. However, for the latter result, more research is needed to make it directly practically applicable. Regardless, our approach already improves upon the state of the art in RRM-ssRNA modelling and is in principle extendable to other types of protein-nucleic acid interactions.

Predicting nuclear G-quadruplex RNA-binding proteins with roles in transcription and phase separation.

RNA-binding proteins are central for many biological processes and their characterization has demonstrated a broad range of functions as well as a wide spectrum of target structures. RNA G-quadruplexes are important regulatory elements occurring in both coding and non-coding transcripts, yet our knowledge of their structure-based interactions is at present limited. Here, using theoretical predictions and experimental approaches, we show that many chromatin-binding proteins bind to RNA G-quadruplexes, and we classify them based on their RNA G-quadruplex-binding potential. Combining experimental identification of nuclear RNA G-quadruplex-binding proteins with computational approaches, we build a prediction tool that assigns probability score for a nuclear protein to bind RNA G-quadruplexes. We show that predicted G-quadruplex RNA-binding proteins exhibit a high degree of protein disorder and hydrophilicity and suggest involvement in both transcription and phase-separation into membrane-less organelles. Finally, we present the G4-Folded/UNfolded Nuclear Interaction Explorer System (G4-FUNNIES) for estimating RNA G4-binding propensities at http://service.tartaglialab.com/new_submission/G4FUNNIES .

Immunofluorescence microscopy of G-quadruplexes and R-loops.

A large variety of non-B secondary structures can be formed between DNA and RNA. In this chapter, we focus on G-quadruplexes (G4) and R-loops, which can have a close structural interplay. In recent years, increasing evidence pointed to the fact that they can strongly influence each other in vivo, both having physiological and pathological roles in normal and cancer cells. Here, we detail specific and accurate methods for purification of BG4 and S9.6 antibodies, and their subsequent use in immunofluorescence microscopy, enabling single-cell analysis of extent and localization of G4s and R-loops.

Identification of bona fide RNA G-quadruplex binding proteins.

RNAs often accomplish their diverse functions through direct interactions with RNA-binding proteins (RBPs) in a sequence- and/or structure-dependent manner. RNA G-quadruplexes (rG4s) are unique secondary structures formed by guanine-rich RNA sequences which impact RNA function independently and in combination with RBPs. Efforts from several labs have identified dozens of rG4 specific RBPs (rG4BPs), although the research is still in the growing phase. Here we present methods for the systematic identification of rG4BPs using a pull-down approach that takes advantage of the chemical modification of guanine bases. This allows abolishing the rG4 structures while still maintaining the base composition intact, and hence helps in recognizing true rG4BPS (in contrast to G-rich motif binders). In combination with other biochemical assays, such an approach can be efficiently used for the identification and characterization of bona fide rG4BPs.

In-gel staining methods of G4 DNA and RNA structures.

G-quadruplexes (G4) are functionally important nucleic acid structures, involved in many cellular pathways. They are often dynamically regulated in cells, which makes detecting them in vivo challenging and dependent on sophisticated technical equipment. Therefore, in vitro studies are commonly performed as a first step to confirm a candidate sequence folds into a G4. Several methods have been developed, each with its individual pros and cons. A highly accessible and quick approach, without the need for specialized equipment, is the detection of G4s in native gels using light-up probes. These molecules become fluorescent after specifically binding to G4s. Several different classes have been discovered, emitting light in various colors, and some possess specificity for certain G4 topologies, which makes them highly versatile tools for G4 visualization. Here, we will explore the general procedure using the light-up probe NMM on RNA G4s and discuss advantages and limitations of this method.

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields.

Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.